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Abstract 
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We study the transmission through different small systems as a function of the coupling strength 
v to the two attached leads. The leads are identical with only one propagating mode £^ in each of 
them. Besides the conductance G, we calculate the phase rigidity p of the scattering wave function 
^>f, in the interior of the system. Most interesting results are obtained in the regime of strongly 
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overlapping resonance states where the crossover from staying to traveling modes takes place. The 
crossover is characterized by collective effects. Here, the conductance is plateau-like enhanced in 
some energy regions of finite length while corridors with zero transmission (total reflection) appear 
in other energy regions. This transmission picture depends only weakly on the spectrum of the 

O 

closed system. It is caused by the alignment of some resonance states of the system with the 
propagating modes £q in the leads. The alignment of resonance states takes place stepwise by 

> 

resonance trapping, i.e. it is accompanied by the decoupling of other resonance states from the 
continuum of propagating modes. This process is quantitatively described by the phase rigidity 
p of the scattering wave function. Averaged over energy in the considered energy window, (G) is 

o 

correlated with 1 — (p). In the regime of strong coupling, only two short-lived resonance states 

O 

survive each aligned with one of the channel wave functions £q. They may be identified with 
^ traveling modes through the system. The remaining M — 2 trapped narrow resonance states are 

^ well separated from one another. 



1 



I. INTRODUCTION 



Resonance phenomena in quantum systems at very low and at very high level density are 
well understood. In the first case, isolated resonances can be seen in the cross section whose 
decay widths are smaller than the distance between them. They are described well by means 
of the discrete states of the closed system. In the second case narrow resonances appear, 
superposed on a smooth background. Here, the resonance states are usually described on the 
basis of statistical assumptions while the background is represented by the so-called optical 
S matrix [T]. Much less understood is the transition between these two borderline cases. In 
this crossover regime, the resonances overlap, and a full (standard) statistical description is, 
generally, not justified. For the results of numerical studies for a realistic system see, e.g., 
Ref . [2] . Characteristic features of the scenario with overlapping resonances are the avoided 
crossings of the resonance states in the complex plane which may change considerably the 
spectroscopic properties (such as positions and widths) of the states [3]. In order to describe 
these changes under realistic conditions (especially in the few-channel case), both the real 
and imaginary parts of the coupling via the continuum have to be taken into account. In the 
standard statistical approach, only the real part is considered. In the recent paper only 
the imaginary part of the coupling via the continuum is taken into account. Here, serious 
deviations from the standard statistical approach are found. Strong deviations are found 
also when both the real and imaginary parts are considered [2]. Generally, the phases of the 
wave functions are not rigid in the regime of overlapping resonances as shown analytically 
in Ref. [5]. 

Mathematically, an open quantum system can be described best by a non-hermitian 
Hamilton operator that accounts not only for the energies but also for the widths (lifetimes) 
of the resonance states. Spectroscopic studies on the basis of a realistic non-hermitian 
Hamilton operator (with complex coupling term) are performed up to now either by using the 
method of complex scaling [H] or the method of the Feshbach projection operator technique 
[3]. The last method is applied also to the transmission through small devices. The theory 
for this case, presented in Ref. [7], is equivalent to the tight-binding lattice model [Bj. In 
the numerical results, long-range correlations between the individual resonance states have 
been observed. 

Most theoretical studies with non-hermitian Hamilton operators are based on random 
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matrices. Usually, the Pandey-Mehta Hamilton operator H(a) = Ho + aHi is used where Hq 
and Hi are real and complex random hermitian matrices, respectively, and a is a crossover 
parameter [9]. Recent theoretical studies have shown that long-range correlations in the 
eigenfunctions of H(a) appear also in these studies [TU] . The theoretical results are in good 
agreement with experimental ones obtained in an open microwave billiard [TTj . 

In order to describe the intensity fluctuations generated by a monochromatic source in a 
disordered cavity coupled to the environment, the idea of standing waves in almost closed 
systems and traveling waves in open systems is considered in Ref. [12J. In this paper, the 
two borderline cases are related to one another by an interpolation procedure. 

Another characteristic feature of overlapping resonances is the fluctuation picture of the 
cross section. The autocorrelation function shows Ericson fluctuations [13J studied experi- 
mentally first in nuclear reaction cross sections [H]. Recently they are studied experimen- 
tally with high resolution in the photoionization of rubidium atoms in crossed magnetic 
and electric fields [JS] . The experimental results obtained in this study proved not only the 
universality of the Ericson fluctuations. They showed an unexpected high reproducibility of 
the results. This experimental result does not coincide with the usual understanding of the 
nature of Ericson fluctuations. 

In recent theoretical studies on Ericson fluctuations, the importance of interferences be- 
tween different resonance states is discussed. In Ref. (16] it is presumed that Ericson 
fluctuations do not result from interferences, but may be considered as a collective coherent 
resonance phenomenon. In the literature, this collective phenomenon is called resonance 
trapping. It is studied theoretically as well as experimentally p2] in many papers, see 
the reviews [3]. It is a collective effect caused by the decoupling of most of the resonance 
states from the continuum at large coupling strength between system and environment while 
only a few states become short-lived. This interpretation of Ericson fluctuations is called 
misinterpretation in Ref. [TS] due to the result obtained by the author, that the fluctua- 
tions originate from interferences between random-partial-width amplitudes in the regime 
of strongly overlapping resonances. 

It is the aim of the present paper to study the transmission through microwave cavities 
as well as through toy-model systems with a small number of sites (corresponding to a 
small number of levels) in the two-channel case (one in each of the attached leads) without 
using any statistical assumptions. The mathematical basis of our study is the Feshbach 
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projection operator [T] formalism, i.e. the non-hermitian effective Hamilton operator H e s 
(see Sect. II) and its relation to the Green function G c g = (E — H^) -1 . The numerical 
results for microwave cavities (Sect. Ill) as well as for one-dimensional and two-dimensional 
toy-model systems with different positions of the attached wires (Sect IV) are obtained 
by using the tight-binding lattice Green function method [8] providing accurate results for 
sufficiently large numerical grids. We turn our attention especially to the appearance of 
coherent collective effects and long-range correlations in the crossover regime from standing 
to traveling waves, i.e. in the regime of overlapping resonances. Furthermore (Sect. V), we 
show the appearance of Ericson-like fluctuations. We are interested in the structure of the 
autocorrelation function for different values of the displacement and for different degrees of 
resonance overlapping. We comment also on the reproducibility of the results by changing 
the geometry of the system. The results are discussed and summarized in the last section. 

II. SPECTROSCOPIC VALUES OF AN OPEN QUANTUM BILLIARD 
A. Non-hermitian Hamilton operator of the open system 

In an exact description of resonance phenomena by using the Feshbach projection operator 
technique pQ, the non-hermitian effective Hamilton operator 

H cG = H B + £ V BC 1 V CB (1) 

C=L,R h ~ h C 

appears where H B is the hermitian Hamilton operator of the corresponding closed system 
with M discrete states, V B c is the symmetrical coupling operator between system and en- 
vironment, He is the hermitian Hamilton operator describing the propagating modes in 
the environment and C stands for the continuum of propagating modes. In the case of a 
quantum cavity, C consists of the waves in the right (R) and left (L) leads attached to the 
cavity [7] . In this representation, the resonance picture is determined by the eigenvalues z\ 
and eigenfunctions 0a of H cG , Eq. ([I]), which differ, generally, from those of H B not only by 
Ihi(za) 7^ 0. Also Ke(z\) differs from the (real) eigenvalue of H B by the contribution arising 
from the principal value integral of the second term of H eG . 

The M eigenvalues z\ of H eG are complex and energy dependent. They provide the 
positions E\ and widths T\ of the resonance states by solving the fixed-point equations 
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E\ = Re(z\)\ E=E ^ and denning I\ = —2 Im(z\)\ E=E . These values correspond to the poles 
of the S matrix. However, for the S matrix describing a physical scattering process at the 
real energy E of the system (without any poles), the solutions of the fixed-point equations 
are relevant only for narrow resonances. Generally, the eigenvalues z\ with their full energy 
dependence are involved in the S matrix [3]. 

Also the eigenf unctions (fix of H e g are complex and energy dependent. The (fix are 
biorthogonal. Due to the symmetry of H efi , it holds [3], 

(<P*x\<fix>) = 5 x ,x> ■ (2) 

As a consequence, {(fix \ (fix) = A x > 1 and A x — > oo with approaching the branch point at 
which two eigenvalues coalesce, zx = zy [191 EOj I2TJ . 

Considering the coupling strength v as control parameter, the phenomenon of resonance 
trapping [3j [22] appears, i.e. the wave functions of a few resonance states align each with 
one of the channel wave functions while the other resonance states decouple from the contin- 
uum (become trapped). The resonance trapping occurs hierarchically [231 [21] • Some years 
ago, the phenomenon of resonance trapping is studied theoretically [22] and even proven 
experimentally [17] on microwave cavities. 

The eigenfunctions (fix of H e ^ can be represented in the set {(fix} of eigenfunctions of the 
Hamilton operator Hb of the corresponding closed system, 

M 

0A = £ rf AA'0£ (3) 
A' 

with complex coefficients dxx 1 normalized according to \dxw\ 2 1 Sa"=i Ma'A"| 2 - 



B. Scattering wave function and transmission 

The eigenfunctions (fix of H e s are part of the total scattering wave function ^/q that is 
solution of the Schrodinger equation (H — E)^/q = in the total function space consisting of 
the discrete states at the energies E^ of the closed quantum billiard and of the propagating 
modes £^ in the attached leads. Due to the coupling between the two subspaces, the discrete 
states turn over into resonance states at the energies Ex of the open billiard. H is hermitian. 
The m% read [3] 

*g = $ + £n?^f^ (4) 

\ E- z x 



where fl x = [1 + (E + — Hc)~ 1 V] <p\ is the wave function of the resonance state A. Accord- 
ing to Q, the eigenfunctions <p\ of the effective Hamilton operator H e g give an essential 
contribution to the scattering wave function in the interior of the cavity, — *> ^/q, 
with 

^c = 1^ °\e <P\ , c XE = — . (5) 

A E-Z X 

The coupling coefficients ((fi*x\V\£.c) are, generally, complex and energy dependent due to 
the unitarity of the S matrix [27]. This energy dependence becomes important only for 
overlapping resonances. Here one of the resonances may appear as some "background" onto 
which a neighboring resonance is superposed [28] . 
In Refs. [10, 26], the phase rigidity \p\ 2 with 

= Jdr fr(r) 2 = c2ig fdr ([Re^(r)] 2 - [M(r)] 2 ) 
P Jdr |^(r)| 2 6 Jdr ([Re^(r)] 2 + [Im^(r)] 2 ) 

defined at the energy E, is introduced in order to describe phenomenologically the influence 
of some impurity onto the phase of the scattering wave function We use this expression 
in order to get information on the phases of the in the regime of resonance overlapping. 
The amplitude of the transmission through a quantum dot is [7J 

a E-z x 

The eigenvalues zx and eigenfunctions <px of H c r are involved in ([7]) with their full energy 
dependence [Bj. According to (|7]), the transmission is resonant in relation to the real part of 
the eigenvalues of H e ^. The transmission amplitude ^ can be rewritten by means of the 
scattering wave function ([5]), 

t = -2m ($|V|fcg> (8) 

with tyf, being complex, in general. The advantage of this representation consists in the fact 
that it does not suggest the existence of Breit-Wigner peaks in the transmission probability. 
Quite the contrary, the transmission is determined by the degree of alignment of the wave 
function ty® with the propagating modes £q in the leads, i.e. by the value (£c\V\^c)- This 
alignment is important in the regime of overlapping resonances where the identification of 
single resonances and resonance states ceases to be meaningful. The expressions ([7]) and Q 
are fully equivalent. 



C. Correlation measures 



The redistribution processes taking place inside the system under the influence of the 
coupling to the continuum can be characterized by the number N% of principal components 
of the eigenfunctions 4>\. Using ^ this number reads [29] 



1 . ,,12 \d 
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TUP _ t ■ \h I 2 — I^ AA/ ' (Q\ 

A MV M \h 14 ' 1 AA 1 sr M U 12 ' ^ > 

M 1^ A'=l I °AA' I 1^ A" =1 I «A' A" | 

In the limit of equal mixing of the state A with all states A', one gets b\\i = l/y/M for all 
A' and N% = 1. In the opposite case, we have b\\i = 5x,y and iV| = 1/M. In the first case 
the collectivity of the state A is large [U [29] . Almost all components of the wave function 
have the same sign because the wave function of the state A is aligned to the scattering wave 
function £q of the channel C at the cost of the other resonance states trapped by it. In 
the second case the wave functions A are not at all aligned to the and the signs of their 
different components are not the same. For details of the phenomenon of alignment see Ref. 
0. Generally, 1/M < N{° < 1. 

Due to the properties of the scattering wave function ^ are determined by the 
properties of all the wave functions 0a- Using the representation we can write = 

Ea(0aI*c) <Px and defme 



Ne £f =1 i<0?i^>r (10) 

The normalization has to be done separately at every energy E and coupling strength v. 
The value Ne is a certain analogon to (|9]), but for the representation of the scattering wave 
function in the set of eigenfunctions 0f of the Hamilton operator Hb- Usually, Ne is 
called inverse participation ratio. Further, we define the analogous expression 

N x = 1^— (11) 

with = (x\i?%) 0f . 

Another measure for coherence is the autocorrelation function of the fluctuation 5a = 
a — a of the cross section a, 

C{e) = ^-r f 5a(E + e)5a(E)dE . (12) 

(7 \I | J I 

Here, / is the same energy interval used to determine a and 5a = a — a. The Lorentzian 
shape of C(e) for small e is characteristic of Ericson fluctuations [13]. In the overlapping 



regime, the width of the Lorentzian is large compared to the average spacing of the resonance 
states. Such a result has been found in many experimental data, above all in the analysis of 
nuclear reaction cross sections, e.g. [2]. It has been found recently also in high-resolution 
experimental data on the photoionization of atoms [15]. In the present paper, we will analyse 
the fluctuation picture of the transmission and denote the corresponding autocorrelation 
function by Cg(s)- 

III. NUMERICAL RESULTS: TRANSMISSION THROUGH QUANTUM BIL- 
LIARDS, ENERGY-AVERAGED PHASE RIGIDITY AND COLLECTIVITY 

The numerical calculations are performed for chaotic quantum billiards of different type 
in the tight-binding lattice model, for details see Datta [5] and [7]. The phase rigidity p is 
calculated in the representation The coupling strength v between cavity and attached 
leads is the same for both leads. It is varied by means of v = k/kq where kq = 1 is the 
hopping matrix element inside the cavity as well as inside the lead while k is the hopping 
matrix element between cavity and lead. The cavities are small and the number of channels 
is one in each of the two attached identical leads. We trace transmission and phase rigidity 
(averaged over the whole energy window between 11 and 34, in units of the widths of the 
leads) as a function of the coupling strength v between cavity and attached leads. 

In Fig. [Tj we present the numerical results for the transmission (top) through a Sinai 
billiard and for the phase rigidity (bottom) as a function of energy E and coupling strength 
v. The shape of the billiard is shown in Figs. [5]^a,b). Around v « 1, the transmission is 
enhanced and the phase rigidity is reduced at all energies. Analogous results (Fig. [2]) are 
obtained for the transmission through a Sinai billiard with another geometry [Figs. [5](c,d)] . 
Comparing the two figures [T] and [2j we see that the transmission spectrum as well as the 
phase rigidity do change only a little when the geometry of the cavity is varied. This result 
coincides qualitatively with the high reproducibility of the experimental data observed in 
the fluctuations of the spectrum of the atom 85 Rb in strong crossed magnetic and electric 
fields [IS]. 

In Fig. |3j we show the numerical results for the transmission through a billiard of 
Bunimovich type and the phase rigidity as a function of energy and coupling strength. 
In the calculations, the leads are attached in such a manner, see Figs. [5te,f), that the 
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FIG. 1: Conductance G = \t\ 2 (top) and phase rigidity \p\ 2 (bottom) for a Sinai billiard as a function 
of energy E and coupling strength v between billiard and attached leads [30]. The transmission 
and the phase rigidity vary between 1 (white) and (black) . The shape of the billiard is shown in 
Figs. [5ja,b): size x = 3, y = 7, radius R = 1.5, in units of the width of the leads. The calculations 
are performed in the tight-binding lattice model [8 . 

transmission via short-lived whispering gallery modes is supported by the geometry of the 
billiard in the regime of overlapping resonances. The whispering gallery modes are especially 
stable as has been shown in earlier studies [31] and as it can be seen directly by comparing 
Figs. [1] and [2] with [3} 

The correlation between the conductance (G) = (\t\ 2 ) and 1 — (\p\ 2 ) can be seen more 
directly from Fig. El where (G) and (|p| 2 ), averaged over energy, are shown as a function of 
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FIG. 2: The same as Fig. [I] but x = 4, y = 5. The shape of the billiard is shown in Figs. [5]^c,d). 

v for the three cavities whose geometries are shown in Figs. [5^a-f). In all cases, (G) and 
1 — (|p| 2 ) are strongly correlated. Around v = 1 (G) and 1 — (|p| 2 ) are maximum. Here also 



Ne, defined in (10), is maximum. When the parameters of the Sinai billiard are changed 
from x = 3,y = 7, Fig. ma), to x = 4, y = 5, Fig. Qb), (G) and (|p| 2 ) do change only a 
little. 

The situation is however another one for transmission through a Bunimovich cavity in 
the regime of overlapping resonances. When the leads are attached in such a manner that 
the transmission via whispering gallery modes is supported then the transmission is more 
enhanced and the phase rigidity is more reduced [Fig. |4][c)] than in the case of the two 
Sinai billiards considered above. When the leads are attached in another manner to the 
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FIG. 3: The same as Fig. [I] but for a billiard of Bunimovich type to which the leads are attached 
in such a manner that transmission via whispering gallery modes is supported. The shape of the 
billiard is shown in Figs. [5je,f): Radius R = 3 and distance D = 2 between the centers, in units 
of the width of the leads. 

Bunimovich cavity so that they do not support transmission through whispering gallery 
modes, the transmission (averaged over energy) may be very small and, correspondingly, 
the phase rigidity (averaged over energy) remains large. For the geometry shown in Figs. 
|5^g,h), it is (G) < 0.1 and (|p| 2 ) > 0.9. In this case, the whispering gallery modes are 
related, above all, to the reflection. 

In Fig. [5j we show the wave functions <p\ of a few short-lived and long-lived states in 
the four cavities considered. The corresponding eigenvalues z\ calculated at the energy E 
are given m Table |T} Although the wave functions are delocalized in space in all cases in the 
corresponding closed system, the short-lived states might be localized as can be seen from 
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FIG. 4: (color online) Conductance (G) = {\t\ 2 ) (blue full line), energy-averaged phase rigidity 
(\p\ 2 ) (red dashed line) and inverse participation ratio Ne of the scattering wave function 
(green dash-dotted line), averaged over the whole energy window (11 < E < 34), as a function of 
the coupling strength v, [30 . (a) Sinai billiard with x = 3, y = 7, R = 1.5 in units of the width of 
the leads [Figs. |5](a,b)]; (b) the same as (a) but x = 4, y = 5 [Figs. [5|c,d)]; (c) Bunimovich billiard 
with transmission through whispering gallery modes, R = 3, D = 2 in units of the width of the 
leads [Figs. [5^e,f)]. The calculations are performed in the tight-binding lattice model [8]. 

Figs. |5](a), (c), (e) and (g). Only the long-lived resonance states are delocalized inside the 
cavity, see Figs. [5](b), (d), (f) and (h). 
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FIG. 5: (color online) Wave functions <j)\ of some short-lived and long-lived resonance states in the 
Sinai billiard (a,b) with x = 3, y = 7; (c,d) with x = 4, y = 5, and in the Bunimovich billiard (e,f ) 
with leads attached for support of transmission via whispering gallery modes and (g,h) without 
support of whispering gallery modes. The corresponding eigenvalues z\ are given in Table [TJ The 
short-lived states (a,c,e,g) are localized in space and de-localized in energy while the long-lived 
states (b,d,f,h) are de-localized in space and localized in energy. 
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TABLE I: The eigenvalues z\ corresponding to the eigenfunctions <fi\ shown in Fig. [5j calculated 
at the energy E and coupling strength v = 1. 



Figure 


E 




z\ 




5 


a 


25 


20.1228 


- 3.1657 i 




5 


c 


25 


18.4570 


- 2.3354 i 




5 


e 


25 


24.8227 


- 2.4076 i 




5 


g 


20 


22.7141 


- 1.1801 % 




5 


b 


25 


20.6891 


- 0.0617 i 




5 


d 


25 


22.6227 


- 0.0635 i 




5 


f 


25 


32.2046 


- 0.0276 % 




5 


h 


20 


18.3281 


- 0.0713 % 



IV. PHASE RIGIDITY AND RESONANCE TRAPPING IN THE TRANSMIS- 
SION THROUGH FEW-SITES STRUCTURES 

The results shown in Fig. [4] show some collective features involved in the scattering wave 
function \1>^. In order to understand these features better, we provide numerical results for 
simple structures with a small number of sites. First, we consider the one-dimensional case 
(chain) described by the tight-binding Hamilton operator 

= -£*ib'>0' + i| + cc (13) 

where j runs over all sites of the system, i.e. over those of the left and right leads as well 
as over those of the one-dimensional box consisting of N sites. The tj are the hopping 
matrix elements. Details of the model are given in [Tj. The box is opened by varying the 
hopping matrix elements between leads and box in the following manner. We define tj = v l 
if j = j in — 1; tj = vr if j = j out ; tj = 1 in the leads and tj = t otherwise. In the case 
t = 1, the eigenvalue spectrum E n = — 2t cos(j^ L ) of H B is distributed over almost the 
whole propagation band E(k) = — 2 cos(k). Our calculations are performed for symmetrical 
cases with — Vr = v. The phase rigidity is calculated from ([6]). 

In Fig. |6j we show results of calculations obtained for a chain of length M = 10 with 
to = 1. In the sub-figure |6](a,b) the conductance G = \t\ 2 and the phase rigidity \p\ 2 are 
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FIG. 6: (a): The conductance G = \t\ 2 of a one-dimensional 10-sites chain, to = 1, with wires 
attached at the ends of the chain (j- m = 1, j out = 10) as a function of energy E and coupling 
strength v between the chain and the two wires [30]. The positions of the long-lived states at 
large coupling strength v differ from the positions E? (denoted by white circles), (b): The phase 
rigidity \p\ 2 for the same 10-sites chain, (c): Transmission (G) (full line) and y/1 — (\p\ 2 ) (dotted 
line), and (d): N x (full line) and Ne (dotted line), as a function of the coupling strength v for the 
10-sites chain. The values (G), {\p\ 2 ), N x and Ne are averaged over energy in the whole energy 
window —2 < E < 2. (e) E\ and (f) T\ as a function of the coupling strength v. The E\ and T\ 
are calculated at E = 0.5. All calculations are performed in the tight-binding model [7J. 

shown as a function of energy E and coupling strength v. These figures are very similar 
to Fig. [3] for the transmission through a quantum billiard with convex boundary in which 
the transmission occurs mainly via whispering gallery modes. In the one- dimensional chain 
(Fig. [6]), we see threshold effects at E = — 2 and E = 2 where the widths of the resonance 
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states approach zero. These threshold effects are excluded in Fig. [3] where the energy range 
13 < E < 34 is considered and the thresholds are at about E — 10 and 40. We conclude 
therefore that the one-dimensional chain mimics the transmission through the Bunimovich 
cavity occurring, at v = 1, mainly via whispering gallery modes. The transmission has a 
plateau at v cv ~ 1. In the case of the Bunimovich cavity, dips appear due to the long- 
lived trapped resonance states. Also the phase rigidity |p| 2 of the chain [Fig. |6](b)] shows a 
behavior similar to that obtained for the Bunimovich cavity. In the case of the chain, the 
transmission occurs clearly via traveling modes. 

For further discussion, we show additionally in Fig. [6] some characteristic values discussed 
in the foregoing sections. Also in the one-dimensional chain, we see the anti-correlation 
between conductance G = \t\ 2 and energy-averaged phase rigidity \p\ 2 [Fig. |6](c)] that is 
observed in the billiards of Sinai and Bunimovich type (Figs. [4]). The values G = \t\ 2 and 
yT — |p| 2 , averaged over energy, are almost the same for all values v. Furthermore, Ne is 
enhanced in the critical region around v ~ I [Fig. |6^d)] in a similar manner as in Fig. [4] for 
the cavities. Additionally, we show N x as a function v. It is also enhanced in the critical 
region. 

In the chain, N x is smaller at v ^> 1 than at v <C 1 according to the fact that the wave 
functions of the long-lived states appearing at large v are, as a rule, more chaotic than the 
original ones at small v (see e.g. [32]), i.e. they are distributed in x with more or less equal 
probability over the whole cavity in contrast to the original states at small v. As to Ne, it 
shows the opposite behavior, i.e. it is larger at v ^> 1 than at v -C 1. The reason is that 
the short-lived states are strongly correlated, i.e. the number N% of principal components 
defined in (|9]), is large for every short-lived state [31 [29]. Due to this fact, Ne remains 
comparably large when short-lived states are formed. 

In the lower part of Fig. [6j we show the energies E\ and widths T\ of the 10 resonance 
states as a function of the coupling strength v. We see the formation of two short-lived states 
in the critical region v ~ 1 according to the fact that two (identical) wires are attached. The 
widths of all states are maximal in the center of the energy window due to the condition 
T\ —>■ for all A in approaching the thresholds E — > ±2. The short-lived states arise 
therefore from the middle of the spectrum. Their contribution is shifted to energies out of 
the window —2 < E < 2 when E ^ 0, see Fig. |6^e) calculated at E = 0.5. We see therefore 
only 8 resonances in the transmission when v > 1. Collective effects are maximal around 
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the critical value v CI = 1 being somewhat smaller than the value v = 1.2 at which the widths 
bifurcate. At v CT , many resonance states are almost aligned with the propagating modes £q 
in the leads, while at v = 1.2 two resonance states are aligned at maximum. 

In order to get a better understanding for the appearance of collective effects in the 
critical region around v CI we show in Figs. |7| to 10 the results of calculations performed 
for structures being more complicated than a chain. We performed the calculations for 
different positions j- m and j out at which the wires are attached, and for different numbers M 
of resonance states. We choose t = 0.05 (instead of t = 1 in Fig. [6]). Due to the small 
value of to, the resonance states lie in the narrow energy range —0.1 < E < 0.1, threshold 
effects are excluded and the eigenvalues z\ are almost independent of energy. 



In the sub-figures (a) of Figs. [7| to 10, we show the transmission G = \t\ 2 against 
coupling strength v and energy E for different one-dimensional systems. The areas with 
0.98 < \t\ 2 < 1 are encircled by red full lines and those with < |t| 2 < 0.02 by white dashed 
lines. In these areas, the transmission at a fixed value v shows a plateau of about 1 and 0, 
respectively. While the positions in energy of the plateaus \t\ 2 « 1 depend on v, the energies 
of the plateaus \t\ 2 ~ do almost not vary with v. They are corridor-like and divide the 
resonances into different groups. The number of corridors is equal to j- m — 1 corresponding to 
an arrangement of the M resonances in j[ n groups. In each group, plateaus \t\ 2 ~ 1 appear 
in the critical region of v, while the resonances cause transmission peaks of Breit-Wigner 
shape at v <C v cr and v ^> v CT . The number of transmission peaks is equal to the number M 
of resonance states when w < « CI , but M — 2 when v 3> v cr due to the fact that two channels 
are open with which two of the resonance states align eventually in the critical region. 

In Figs. to |H the wires are attached symmetrically each to one of the internal sites, 
jout = M — j in + 1, j in > 1. In difference to the case represented in Fig. |6j the incoming 
wave is split into two parts: both parts move in different directions each to one of the 
borders of the chain where it will be reflected. Then the two parts interfere and give rise to 
a complicated transmission picture. As a consequence of the interferences, the transmission 
picture does almost not depend on the number of the sites (M = 10 in Figs. [7] and |8j 
but M — 11 in Fig. [9J. More important than the total number of sites is the distance 
of the attached wires from the border of the chain. When the distance is an odd number 
(M- m = 2, M out = M — 1 in Fig. [7]) there is a zero transmission corridor around E — 0, 
while the transmission is maximal around this energy when the distance is an even number 
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(M in = 3, M out = M — 2 in Figs. [8]and[9]). The transmission picture shows more structures 
when the symmetry is violated (Fig. |l0~|for M — 11 with M in = 2, M out = M - 2). 

The regions of extremely low conductance (less than 0.02, corresponding to transmission 
gaps) and enhanced transmission (more than 0.98, corresponding to perfect conductance 
bands), marked in Figs. [7] to [9] by white dashed lines and red solid lines, respectively, are 
extended over comparably large energy regions. A similar result ("transmission gaps and 
sharp resonant states") has been found in recent numerical calculations for the transport 
through a simple mesoscopic device [S3]. We underline that both the transmission plateaus 
with maximal transmission and the corridors of transmission zeros are the result of redis- 
tribution, interference and alignment processes. They can not be described by isolated 
Breit-Wigner resonances or isolated bound states in the continuum. 

The redistribution processes taking place in the critical region around v CI = ^/to of the 
coupling strength v can be traced by comparing the sub-figures (b) of Figs. [7] to [To] with the 
corresponding sub-figures (c). In the sub-figures (b), the coupling coefficients \v(ri)\ of the M 
discrete states [eigenstates of H B , Eq. (fll)] to the propagating modes £^ are shown. They are 
different for the different positions j m of attached wires and different numbers M of resonance 
states. Due to the symmetries involved in the considered chains, differences appear between 
the cases with an even and an odd number M of resonance states. In the first case [M = 10, 
Figs. [7] and [8], the coupling coefficients v(n) (1 < n < M) of all M states to the continuum 
of propagating modes are different from zero while there is at least one state with v(n) = in 



the last case [M = 11, Figs. [9] and 10 . Nevertheless, the eigenvalue trajectories [sub-figures 
(c)] are similar for the different cases with different M. They show the typical resonance 
trapping behavior known from studies on many different systems, see Ref. [3j. In the 
regime of overlapping resonances and two channels, the redistribution of the spectroscopic 
properties of the system takes place due to the alignment of two resonance states each with 
one of the two propagating modes (input and output channels). The redistribution occurs 
hierarchically [231 El]: a resonance state will be trapped by overlapping with a neighboring 
state when its width is somewhat smaller than that of the neighboring state. The widths 
of trapped states cease from increasing with increasing coupling strength v. Instead, they 
mostly decrease and the positions in energy of the trapped states do almost not change with 
further increasing v, for details see Ref. [3J. 

The formation of zero-transmission corridors is directly related to the existence of trapped 
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resonance states. They appear independently of the existence of states whose coupling 
strength v(n) is zero. Between the corridors, the redistribution processes create plateaus 
with an enhanced transmission. These plateaus are related directly to the avoided level 
crossings appearing in the regime of overlapping resonances. Unlike the zero-transmission 
corridors, the positions of the plateaus with maximal transmission depend on the coupling 
strength v. 

In the sub-figures (d), the energy-averaged transmission (G) = (\t\ 2 ) and the phase 
rigidity (|p| 2 ) of the total scattering wave function are shown. Like in the calculations 
for the microwave cavities, the transmission is enhanced and the phase rigidity is reduced 
in the critical region where the spectroscopic redistribution takes place. 

The calculations for two-dimensional systems are performed in an analogous manner as 
those for one-dimensional systems, for details see Ref. [7j. The results are similar to those 
obtained for the one-dimensional systems. The main difference between the two cases is that 
there are more possibilities to influence the transmission picture in the two-dimensional 



systems than in the one- dimensional systems. An example is shown in Fig. 11 It is 
interesting to see that the areas with enhanced transmission may be shifted and that the 
corridors with zero-transmission zeros may be relatively broad. 



Summarizing the results shown in Figs. [6] to 11 for different toy models, we state the 
following. The enhancement of the transmission seen in realistic systems in the crossover 
from the weak-coupling to the strong-coupling regime (Sect. Ill) is caused by spectroscopic 
redistribution processes taking place in the system under the influence of the coupling to 
the continuum. These redistribution processes are correlated with a reduced value of the 
phase rigidity and may be traced back to the well-known resonance trapping phenomenon 
occuring in open quantum systems at high level density when the resonance states overlap. 

V. FLUCTUATION OF THE TRANSMISSION PROBABILITY 

We performed an analysis of the transmission through the considered cavities by means 
of the autocorrelation function of the transmission fluctuations 5G = 5\t\ 2 = \t\ 2 — (\t\ 2 ), 



see Eq. (12). In the calculations, / = 23 is the energy interval used to determine (\t\ 2 ), 
i.e. 11 < E < 34 as shown in Figs. [TJ [2] and [3j The area of all three cavities is of the 
same order of magnitude. The number of states is about 35 in the considered energy region 
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11 < E < 34. The average spacing of the states is about 0.65. There is one channel in each 
of the two identical attached leads. 



In Fig. 12, we show the autocorrelation function Cg(s) for the transmission through the 
three cavities (considered in Sect. Ill) in the regime of overlapping resonances with avoided 
level crossings where the transmission is enhanced and the phase rigidity is reduced. At small 
e, we see the Lorentzian-like shape characteristic of Ericson fluctuations |13j . When v = 1, 
the width of the Lorentzian is larger than the average spacing. Comparing the function 
Cq{e) of the two Sinai billiards with that for the transmission through the Bunimovich 
cavity via whispering gallery modes we see that the width of the Lorentzian is larger for 
the Bunivovich billiard than that for the two Sinai billiards. In all three cases, the width 
becomes smaller with decreasing as well as with increasing coupling strength, v = 0.9 and 
1.1, respectively. Enlarging v further, we arrive in the regime of separated time scales. Here, 
Cq{e) falls off in a smaller region e than in the regime with overlapping resonances. That 
means, the correlation between the long-lived states is small at large v and is comparably 
in value with the correlations between almost isolated resonance states (at small v). 



As can be seen from Fig. 12, the function Co(e) is still large at larger e, especially in 
the case of the Bunimovich cavity. This result points to long-range correlations existing 
between the different states, above all between different short-lived states. The oscillations 
differ from one another in the three considered cases. Related to this statement, we point to 
another result of our calculations: in the autocorrelation function Cq{s) at f ^> 1, we have 
small "peaks" (instead of oscillations) that are related to distant narrow resonance states. 

In Fig. [l3j the autocorrelation function Cq versus energy E and coupling strength v 
is shown for the transmission through the 11-sites chains, represented in the insets of the 



sub-figures (c) of Figs. [9] and 10 with, respectively, symmetrically and non-symmetrically 
attached leads, for four different values v around the critical value v cr . Both cases show 
the Lorentzian-like peak with maximum in the critical region v — v CT — ~ 0.2. The 
width of the peak increases with v as long as v < v CI . This result agrees with theory [13J. 
Approaching however v CI , the width of the Lorentzian is not well defined. For v > v CI , the 
mean width of the long-lived resonance states decreases as v~ 2 (see e.g. Fig. [6]), while the 
width of the Lorentzian decreases much less. The structure of the autocorrelation function 
is more pronounced in the non-symmetrical case than in the symmetrical one. Altogether, 



the results agree with those obtained for the Sinai and Bunimovich cavities (Fig. 12). 



20 



The fluctuations in the transmission picture result, of course, from interferences. However, 
the alignment of some of the wave functions with the propagating modes in the wires plays 
an important role. It is expressed by the phase rigidity as shown in Sects. Ill and IV. The 
numerical results on the fluctuations demonstrate the weak dependence of the structure of 
the transmission picture on external parameters such as the shape of the billiard in Sect. Ill 
or the number of sites in the chains in Sect. IV. They are a hint to the existence of collective 
effects, that are large in the regime of overlapping resonances. 

We mention here the results of three completely different experimental studies that qual- 
itatively agree with our theoretical results on the collective effects observed in the conduc- 
tance in the regime of overlapping resonances. First, the directly measured average life time 
of compound nucleus states in 58,60 Ni(p,p') is significantly smaller than that of the observed 
structures in the excitation function. Moreover, it decreases with increasing bombarding 
energy of the projectile and is smaller than the average distance between the states [34]. 
Secondly, a surprisingly high reproducibility of the fluctuation picture has been observed in 
high-resolution experimental results on Ericson fluctuations in atoms [T5]. Thirdly, long- 
range correlations between individual resonance states have been found experimentally in 
an open microwave billiard 

VI. DISCUSSION OF THE RESULTS AND SUMMARY 

In this paper, we studied the transmission through different mesoscopic systems with 
M states as a function of the coupling strength v between system and attached leads and 
wires, respectively. We considered the two-channel case, i.e. one propagating mode in each 
of the two attached identical leads (wires). The numerical calculations are performed by 
using the Green function method (without any statistical assumptions). The results show 
the following features. At small v, the transmission picture consists of M isolated resonance 
peaks (corresponding to standing waves) while M — 2 narrow resonances are superposed by 
a smooth "background" transmission (caused by the two traveling modes) at large v. Most 
interesting are the results in the regime of strongly overlapping resonance states where the 
crossover from standing to traveling waves takes place. Here, the transmission is enhanced. 

According to our results, the crossover occurs by means of the alignment of the wave 
functions of the individual resonance states each with one of the propagating modes ^c=r 
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and £c=L i n the two leads. The alignment of the resonance states occurs stepwise. It is 
accompanied, at each step, by a certain decoupling of at least one neighboring resonance 
state from the continuum of propagating modes (resonance trapping). For the mathematical 
foundation of the alignment phenomenon in open quantum systems see Ref. [5]. 

Resonance trapping in the regime of overlapping resonance states occurs hierarchically as 



can be seen from Figs. [6] to ll_ It is a collective effect to which, eventually, all the states in 
the considered energy region contribute. Therefore, also the reduction of the phase rigidity 
and the enhancement of observable values, such as the transmission, are collective effects in 
the considered energy region. The reduction of the phase rigidity and the enhancement of 
the transmission are maximal when a large number of individual resonance states is almost 
aligned with the propagating modes in the leads. This happens at a coupling strength v 
being somewhat smaller than the critical one at which the number of aligned resonance 
states is exactly equal to the number of propagating modes (channels). This situation is 
illustrated best by the whispering gallery modes appearing in a cavity of Bunimovich type, 
Fig. [3j see also Figs. [6] to 11 



Most interesting result of our numerical studies is that, in the regime of overlapping 
resonances, the conductance (G) is enhanced and correlated with 1 — (|p| 2 ). For certain 
values of the coupling strength v between system and environment, the conductance is 
plateau-like in some finite energy regions AE. Furthermore, corridors with zero transmission 
(total reflection) may appear as a function of the coupling strength v due to destructive 
interferences between neighboring resonance states. The transmission picture with plateaus 
of maximal transmission and corridors of zero transmission is only weakly influenced by the 
spectrum of the closed system. It results mainly from collective interference effects. 

The relation of the results obtained by us to the idea of standing waves in almost closed 
systems and traveling waves in open systems [12] is the following. In an almost closed 
system it is H e g « Hb, while in a strongly opened system the main term of H e s is the 
coupling term J2c=l,r Vbc{E + — Hc)~ x Vcb via the continuum, see Eq. ([I]). The first case 
is realized in the regime of non-overlapping resonances while in the second case a smooth 
background (arising from the short-lived resonances) is superimposed by long-lived narrow 
(non-overlapping) resonances. Furthermore, the monochromatic source considered in Ref. 
[T2] corresponds to the one-channel continuum represented by the propagating modes in 
(OJ). It is therefore not astonishing that the results obtained from ^ for the transmission 
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amplitudes in the two borderline cases fit well to the picture described in Ref. [12] for waves 
propagating in a random medium. However, the crossover between the two borderline cases 
is described differently in the two methods. In [T2] a simple interpolation between the two 
regimes is proposed. In our formulation, however, the crossover between the two regimes is 
calculated in a straightforward manner. The two equivalent exact expressions and @ 
are used for, respectively, the two limiting cases and the crossover regime with overlapping 
resonances. As a result of these calculations, we received the result that the crossover regime 
is dominated by coherent collective phenomena. 

The collective effects in the regime of overlapping resonance states influence also the 
autocorrelation function of the fluctuations of the transmission probability. Long-range 
correlations do appear. The autocorrelation function does not depend sensitively on external 
parameters such as the shape of the Sinai billiard. The high reproducibility of the fluctuation 
picture observed experimentally on atoms [15j fits into this picture. Furthermore, our results 
coincide with those obtained in Ref. [16] : in the few-channel case, Ericson-like fluctuations 
are determined by collective coherent interferences between the resonance states. 

Summarizing we state the following. The transition from the weak-coupling to the strong- 
coupling regime is controlled by a redistribution of the spectroscopic properties of the system 
under the influence of the environment to which the system is coupled (by attaching leads 
to it). The redistribution consists in the alignment of a few resonance states each with 
one of the propagating modes in the leads. The alignment is a collective effect to which 
all resonance states in the considered energy region contribute: it occurs hierarchically by 
trapping neighboring resonance states which are more weakly coupled to the continuum of 
propagating modes than those which align. Eventually, the number of aligned resonance 
states is equal to the number of propagating modes in the leads. Due to their short life time 
(large decay width), the aligned resonance states form the background scattering on which 
the trapped resonance states appear as narrow resonances. Since they are aligned with the 
propagating modes £<f in the leads, these short-lived resonance states may be identified with 
traveling modes. 
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FIG. 7: (color online) The transmission through a one-dimensional chain of 10 sites, see the inset 
of sub-figure (c). The arrows show the sites at which the two (left and right) wires are attached. 
The hopping matrix elements inside the chain are 0.05 while those inside the attached wires are 1. 
(a) Conductance G versus incident energy E and coupling matrix element vl = vr = v, [30]. The 
white dashed lines show regions of extremely low conductance (less than 0.02, corresponding to 
transmission gaps) while the red solid lines show regions with enhanced transmission (more than 
0.98, corresponding to perfect conductance bands). The positions of the bound states are shown 
by white stars, (b) The coupling matrix elements \V n \ = \i/j n (2)\ = \ifj n (9)\ between the wires and 
the chain, ip n (j), are the eigenfunctions of the "closed" chain described by Hb, see Eq. ([!]). (c) 
Evolution of the complex eigenvalues z\ of the effective Hamilton operator, Eq. (IT]), with increasing 
coupling strength v. The are calculated in site representation according to Datta [8]. The values 
at v = are shown by red circles, (d) Evolution of the mean conductance (G) (red broken line) 
and energy-averaged phase rigidity (\p\ 2 ) (blue full line) with v. Note the hierarchical trapping of 
resonance states shown in (c): The resonance states 1 and 5 (10 and 6) are trapped already at 
small v. The state 3 (8) first traps the state 2 (9) and becomes trapped by the state 4 (7) at larger 
v. The two states 4 and 7 align with the two states in the two channels. 



26 





FIG. 8: (color online) The same as Fig. [7]but \V n \ = |V>n(3)| = |VVi(8)| in (b). Note that resonance 
trapping shown in (c) starts at larger values of v than in Fig. [7J Resonance state 2 (9) traps first 
the resonance states 1, 3, 4 (10, 8, 7) and becomes trapped by state 5 (6) at larger v. The two 
states 5 and 6 in the center of the spectrum align with the two states in the two channels. 
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FIG. 9: (color online) The same as Fig. [7]but for a chain of 11 sites and |V(n)| = 1^(3)1 = |"0 n (9| 
in (b). Note that |Im(z^)| of the resonance state 4 (8) does not increase with increasing v even at 
small v, i.e. this state is trapped already at very small v. The state 5 (or 7) traps first the states 
1, 2, 3 (11, 10, 9) and becomes trapped by the state 6 (in the center of the spectrum) at larger v. 
Finally, also in this case |Im(^)| of 2 states is large at large v, i.e. two states are aligned each to 
one of the states in the two leads. 
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FIG. 10: (color online) The same as Fig. [7]but for a chain of 11 sites and |Vl(w)| = |^ n (2)| (blue 
circles), |Vh(n)| = |^ n (9)| (red crosses) in (b). The regions with enhanced transmission (more than 
0.98, corresponding to perfect conductance bands) are small and not shown. Note the hierarchical 
resonance trapping also in this case. The resonance states 1, 4, 5 (11, 8, 7) are trapped at relatively 
small v while the state 2 (10) first traps the neighboring states and then becomes trapped by the 
state 3 (9) at larger v. With further increasing v, the state 3 (9) becomes trapped by the state 6 
in the center of the spectrum. Finally, there are two states (6 and 3 or 9) with large widths the 
wave functions of which are aligned each to one of the two wave functions in the two leads. 
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FIG. 11: (color online) The same as Fig. [7] but for a rectangle of 6 sites along the x-direction 
and 3 sites along the y-direction, see the inset in (c). The coupling matrix elements |Vi(m, n)| = 
\ipm,n(i = 2, j = 2)|, \Vji(m, n)\ = \ip m) n(i = 5,j = 2)| coincide because of the symmetry of the 
system, the i^m,n(hj) are the eigenfunctions of the "closed" rectangle described by Hb, see Eq. |T]). 
The regions with enhanced transmission (more than 0.98, corresponding to perfect conductance 
bands) are small and not shown. Note that |Im(z,\)| of the resonance state 4 (10) does not increase 
with v even at very small v. Resonance trapping occurs hierarchically first around the center of 
the spectrum and then at its edges. Finally, two resonance states align to the two scattering states 
in the two leads. 
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FIG. 12: (color online) The autocorrelation function Cq(e) for the transmission through the Sinai 
billiard with x = 3,y = 7 (a), with x = 4, y = 5 (b) and through the Bunimovich billiard via 
whispering gallery modes (c). In all three cases, the autocorrelation function is calculated for the 
whole energy window 11 < E < 34 and shown for v = 0.9 (red dashed line), v = 1 (blue full line) 
and v = 1.1 (green dash-dotted line). 
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FIG. 13: (color online) The autocorrelation function Cg(s) for the conductance G = \t\ 2 of the 
11-sites chain with, respectively, symmetrically (top) and non-symmetrically (bottom) attached 
leads [see the insets of the sub-figures (c) of Figs. [9]and[To] for four different values of the coupling 
strength v: v = v CT — 0.1 (dashed blue curves), v = v cr (solid red curves), v = v cr + 0.1 (dotted 
green curves), and v = v cr + 0.2 (dash-dotted black curves). 
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